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We use extensive first principle simulations to show the major role played by interfaces in the 
mechanism of phase separation observed in semiconductor multifunctional materials. We make an 
analogy with the precipitation sequence observed in over-saturated AlCu alloys, and replace the 
Guinier-Preston zones in this new context. A new class of materials, the a phases, is proposed to 
understand the formation of the coherent precipitates observed in the GeMn system. The interplay 
between formation and interface energies is analyzed for these phases and for the structures usually 
considered in the literature. The existence of the a phases is assessed with both theoretical and 
experimental arguments. 
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INTRODUCTION 

Multifuuctioual materials have drawu a lot of atteu- 
tiou iu receut years because of their intrinsic scientific 
interest and for the huge technological possibilities they 
offer. Indeed these materials simultaneously exhibit two 
or more characteristics that are usually incompatible: 
ferroelectricity and ferromagnetism for multiferroics, un- 
correlated thermal and electrical conductivity for ther- 
moelectrics, or semiconductor (SC) character and fer- 
romagnetism. Some homogeneous materials are avail- 
able, however they usually exhibit weak properties or are 
limited to very low temperatures^. In order to open a 
broader range of opportunities, heterogeneous materials 
have been proposed as an efficient way to couple strong 
responses at room temperature^. Among the multifunc- 
tional materials, those based on semiconductors yield 
great promises, and similar enhancements are expected 
in nanostructured media for thermoelectricity^ and mag- 
netic semiconductivity^. Regarding this last issue, it has 
been shown that depending on the solubility limit and the 
growth temperature, magnetic impurities tend to make 
a solid solution and form dilute magnetic semiconductor 
(DMS), or gather in solute-rich precipitates via a phase 
separation process^. The formation of elongated nano- 
clusters observed in several systems such as Al(Cr)N 
(Ref. 5) or (Ge,Mn)-nanocolumns^ is a way of getting 
the attractive heterogeneous materials discussed above. 
Controlling the growth thanks to a detailed understand- 
ing of the phase separation mechanism, is a prerequisite 
for their efficient use to design new multifunctional ma- 
terials. 

As a new perspective, it is well known from physical 



metallurgy that phase separation into an over-saturated 
solid solution can occur via a three-step mechanism, for 
example in the AlCu system^: i) the solute atoms diffuse 
in the host matrix to form Guinier-Preston (GP) zones, 
which are solute-rich zones (also referred to as 0") that 
keep the host matrix atomic structure; ii) these zones 
transform to a secondary metastable phase {0' ^ isostruc- 
tural to CaF2) that is fully coherent with the host matrix, 
thus minimizing interface energy ; iii) a last transforma- 
tion occurs with the formation of the stable AI2CU phase 
6 which is structurally incompatible with the host matrix. 
It is clear in this sequence that interfaces play a major 
role, especially in the stabilization of the metastable 9' 
phase. 

In this article, we show that this interface- driven mech- 
anism is also relevant in the case of magnetic semicon- 
ductors, but with specificities not taken into account up 
to now. In order to get a deeper insight into the problem, 
we consider the GeMn system as a prototype, in which 
either dilute Mn atoms^'^, Mn-rich nanocolumns^'^^'^^ 
(NC), or precipitates of the stable phase^'^^ are observed 
depending on the growth conditions. 

The AlCu precipitation sequence can be applied in the 
case of the GeMn system with the phase being the 
Mn^Ges compound while the 0^ and 0^^ phases as re- 
vealed in the nanocolumns remain unknown. This anal- 
ysis is more complete than the simple spinodal decompo- 
sition scheme^ that only focuses on the first stage of the 
former sequence leading to a coherent solute-rich solid 
solution (i.e. DMS -^ 0"). Still, the precipitation se- 
quence needs to be slightly adapted. On one hand, in 
the case of AlCu, one deals with compact metallic phases 
with atomic structures and chemical behaviors basically 



all compatible in essence, leading to easy mixing, for- 
mation of solid solutions and natural interfaces. On the 
other hand, in the case of dilute semiconductors, the in- 
volved phases display non compact structures (diamond 
or wurtzite) for the host, and much more compact phases 
for precipitates, with very little structural and chemical 
compatibility to one another ^^, leading to strong chemi- 
cal segregation and incoherent interfaces. 

The thorough study of two types of GeMn materials 
in bulk conditions, followed by the introduction of their 
interfaces with the semiconductor matrix, allows us to 
unravel not only the mechanisms behind this variety of 
precipitates, but also the atomic structure of the GeMn 
nanocolumns. This yields a new understanding and po- 
tential control over these complex and very promising 
heterogeneous systems. 



It can be checked that the formation energies [Eq. (1)] of 
the Ge lattice and the MnnGeg compound (or a mixture 
of them) are zero, since they are the chosen reservoirs 
of atoms in the semi-grand canonical ensemble. Note 
that a less natural, but reasonable, choice would have 
been to choose the compound MusGes as a reference, 
since it has a composition close to that of MnnGeg and 
has a simpler unit cell. The difference on the formation 
energies of the compounds studied here would only have 
been of the order of 10 meV. 

This methodology has proven to well reproduce both 
the properties and the stability of all known (Ge, Mn) 
compounds^'^, despite the necessary and usual approxi- 
mations used. 



II. THE SOLID-SOLUTION HYPOTHESIS 



I. METHOD 

We use the density functional theory (DFT) frame- 
work, within the projector augmented-wave approach 
(PAW) and the generalized gradient approximation 
(GGA) for the exchange-correlation energy, as imple- 
mented in the abinit code^^'^^. The details are those 
used in Ref. 13. In particular, the pseudopotentials have 
been generated with the atompaw code^^. They corre- 
spond to valence electron states 4^ and 4p for Ge and 3d 
and 4s for Mn. The /c-point meshes and the cutoff en- 
ergy ensure an error on the total DFT energy lower than 
0.1 eV/atom. 

Since we want to compare the stability of phases with 
different values of the Mn concentration, xun = 1 — ^Ge, 
we calculate their formation energies per atom: 
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where e is the energy per atom found with our first- 
principle calculations, /iMn and /iQe are the chemical po- 
tential of Mn and Ge, respectively. When manganese 
atoms are introduced in a germanium matrix, the ul- 
timate system obtained after annealing is composed of 
MniiGeg grains embedded in the diamond lattice of ger- 
manium. Therefore, these crystals, MnnGeg and Ge, are 
the two reference systems we have considered to deter- 
mine /iMn and /iQe- When two phases 7 and 5^ with dif- 
ferent compositions, are coexisting, the chemical poten- 
tials can be derived by considering the thermodynamic 
semi-grand canonical potential as described in Ref. 13: 
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When 7 is the Ge diamond lattice and 6 is the MnnGes 
compounds, we get: 



I^Ge = e^ 



MMn = (19e^-8eT)/ll 
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We start with the assumption accepted so far^'i^'^i'i'' 
that the binding energy between magnetic impurities in 
the diamond structure drives the formation of a solid 
solution inside the nanocolumns^^. Those are thus sup- 
posed to be domains where the concentration of Mn im- 
purities is locally increased, yet with little effect on the 
diamond structure. Since the Mn atoms could incorpo- 
rate substitutionally or interstitially into the diamond Ge 
crystal lattice, we test the solid-solution assumption by 
considering a class of materials based on both Mn de- 
fects (see Fig. la). More specifically, we use cubic cells 
of 5.658 A of side (8 atoms of Ge), and consider all the 
possible combinations of substitutional and tetrahedral 
interstitial Mn atoms up to a concentration of 60%. In 
the following, we name nTmS a solid solution with n 
Mn in tetrahedral interstitial sites and m Mn in substi- 
tution per cubic unit cell, hence with a concentration 
^Mn = {n -\- m)/{n -\- 8) where n < 4 and m < A. The 
restriction to these only two defects is based on obser- 
vations by transmission electron microscopy showing a 
clear crystalline continuity between Mn-rich regions and 
diamond Ge along the [001] direction, thus discarding 
non crystal- aligned defects. In addition, our simulations 
of six different point defects involving one Mn atom and 
of nine associated Mn dimers confirm the energetic rele- 
vance of only the substitutional and tetrahedral intersti- 
tial defects (see appendix A). Finally, 27 non-equivalent 
configurations are tested (some examples are shown in 
Fig. lb). Ferromagnetic and antiferromagnetic orders 
have been considered, and the antiferromagnetic config- 
urations turn out to have the lowest energies. 

The evolution of the formation energy Ey with the 
Mn content is reported on Fig. 2. First, we can see a 
steep increase of E^ at low concentration (the slope be- 
ing 1.5 eV/Mn atom), in agreement with the low solu- 
bility of Mn in Ge^^. Then at higher concentration, the 
interactions between defects lower the formation energy 
Ey- Consequently, the curve differs from a straight line, 
and the formation energy reaches a maximum around 
400 meV/atom, before decreasing. The most significant 



Mn in diamond Ge 




a (5.66A) 

Diamond Ge 



Substitutional 
Mn 



T: 

Tetrahedral 

interstitial 

Mn 



IT IS 



3T1S 



3T2S 



4T2S 




g-phase family 




a-Ge2Mn 



c 

(-3.7%) 



Simple-cubic Ge 




a (+2.7%) 



d) 



a-Ge2Mn 




Figure 1. (Color online) The two structures tested: solid solution in the left column, intermetallic in the right column, with 
Mn contents XMn indicated at the bottom, a) The diamond Ge lattice and the two most stable isolated Mn defects (S and T). 
b) Four typical combinations of the S and T defects, c) The simple-cubic Ge lattice and the a Ge2Mn phase, isostructural to a 
FeSi2. d) Four different combinations of Mn interstitials, representative of the a-phase family. For clarity reason, the displayed 
structures are unrelaxed. 



result here is this hump-shaped curve with the maximum 
of the formation energy at 35% of Mn, i.e., within the 
concentration range of 15-40% for which the precipitates 
are observed experimentally^' ^^'^^'^^. The absence of a 
minimum of ^F(^Mn) in this experimental range rules 
out a hypothetical stability of a diamond based coherent 
solid solution. 

Accordingly, our calculations show that neither a spin- 
odal decomposition, nor the analog of the first step of an 
AlCu precipitation sequence can occur in the bulk of the 
matrix. In the latter case it would correspond to the so- 
called Guinier-Preston zones (also referred as 0" phase). 
However, subsurface interstitials have been reported in 
Refs. 8 and 20. We propose that this interstitial-rich 
solid solution at the surface could be in fact the 0" phase, 
forming Guinier-Preston zones that will start the next 
precipitation step with a 6^ phase inside nanocolumns. 
In the following section, we propose a structure for such 
a 0' phase. 



III. a-PHASE STRUCTURES 

Since the previous section has shown that GeMn 
solid-solution precipitates are not compatible with the 
£^F(^Mn) curve, we need an equivalent of the 6^ phase of 
AlCu, which is stabilized thanks to its natural interfaces 
with the host matrix. We propose a new class of ma- 
terials based on the insertion of Mn into a simple-cubic 



crystal of Ge (Fig. Ic). 

Such a structure was proposed to describe MnGe4, 
which is stable at high-pressure^^. The exact struc- 
ture was not determined experimentally, but - from sim- 
ilarity with IIg4Pt - it was suggested to involve Mn 
atoms placed every four interstitial positions, which are 
at the center of the Ge cubes. It exhibits a striking 
structural proximity with diamond Ge, with a mismatch 
of only -2.5% and -1.0% along the a and c directions, 
respectively^^. However, the calculated formation energy 
reaches 180 meV/atom^'^. 

Applying the same approach as for the nTmS com- 
pounds, we have tested 45 different atomic configurations 
of Mn interstitials, using various cell sizes to have a com- 
plete energy diagram of all cubic Ge-based compounds. 
We have discovered a whole family of low energy struc- 
tures, based on the Ge2Mn phase shown in Fig. Ic, which 
is isostructural to a FeSi2. We refer to it as the "a-phase 
family" (see Fig. Id). Such a variety in stoichiometry is 
common in the GeMn system, as illustrated by MnsGe 
alloys^'^. The calculated equilibrium lattice parameters of 
the perfect a Ge2Mn phase are a=6=2.90 A and c=5.45 A 
(Fig. Ic), leading to c/a=0.94. Thus, despite its simple- 
cubic crystal structure, this compound exhibits a lattice 
mismatch with bulk Ge lower than 4% in all 3 directions, 
and a cell volume only 1.6% larger, in agreement with 
observations of Ref. 6. Also, a decrease in the Mn con- 
tent is accompanied by a regular decrease of the cell size, 
as expected from the interstitial character of Mn atoms, 
leading to a cell volume at 20% Mn almost equal to that 
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Figure 2. (Color online) Formation energy as a function of Mn 
content XMn, with diamond Ge and GesMnn being the refer- 
ences for the chemical potentials. Solid solution compounds 
are in blue (circles for low concentration DMSs, triangles for 
n Tins' combinations), the a family in red (squares). Hol- 
low markers denote fully-relaxed structures, i.e., both with 
respect to atomic coordinates and cell parameters. Filled 
markers denote structures with the coordinates relaxed while 
the fixed cell parameters are commensurate with the experi- 
mental Ge lattice. The dashed area (green online) shows the 
XMn-range of the experimentally observed nanocolumns. The 
pastel broad lines are guides for the eyes for the solid-solution 
data on one hand, and the a- family data on the other hand 
(respectively, blue and pink online). 



of bulk Ge, with c/a = 0.96. 

Note that the new metastable Ge2Mn compound C16 
described in Ref. 22 is directly related to the "a-phase 
family" but is distorted by a Jahn- Teller- type effect, 
making it strongly tetragonal and incompatible with the 
diamond structure. 

Our calculations show that the a Ge2Mn compound 
exhibits ferromagnetism, with a local moment on Mn 
atoms reaching 1 /iB, and a polarization P0=-80% and 
P2=-83% (with definitions of PO and P2 as in Ref. 23). 

The formation energy of the a phases displays a broad 
minimum around 60 meV/atom, at Mn concentrations 
ranging from 20% to 33%. Above 33%, Mn interstitial 
atoms must be added in the empty interstitial planes of a 
Ge2Mn, at a formation energy cost of about 1 eV per Mn 
atom added. Below 20%, Mn planes of a Ge2Mn become 
less than half filled and the formation energy increases 
rapidly, until the unstable simple-cubic Ge is left. 

Thus, contrary to the nTmS configurations, the a- 
phases family gives rise to a second minimum of en- 
ergy versus stoichiometry compatible with a segregation 
mechanism. Furthermore, the flatness of the minimum 
in the 20 to 33% range, agrees with the different Mn 
concentrations observed in the nanocolumns by Jamet et 
al. (20%-33%)^'^'', as well as in similar precipitates by 
Bougeard et al (15%)^° and by Li et al (17%-30%)i^ 

Since the nanocolumns are embedded in the Ge matrix, 
the elastic energy has to be considered. We do not calcu- 
late exactly this energy, which depends on the nanocol- 
umn distribution, but an upper bound of it. As shown in 



Fig. 2, we have calculated the formation energies of the 
possible structures either with their own lattices param- 
eters, or constrained to be commensurate with the Ge 
lattice (hollow versus filled markers). The last situation 
corresponds to a state in which all the elastic constraints 
are localized inside the columns. In the reality, the ma- 
trix would also adapt to the stress in order to globally 
minimize this energy. Therefore, the energy difference 
between a phase constrained or not constrained by the 
Ge diamond lattice is an upper bound of the total elastic 
energy. These values turn out to be quite small - about 
10 meV/atom - compared to other energy terms, and 
are thus fairly accurate approximations of the values of 
the elastic energy. More precisely, the formation energy 
of a given compound, including the elastic energy due 
to its embedding in the matrix of germanium, is located 
between the hollow and filled markers in Fig. 2. 



IV. ROLE OF INTERFACES 

The structural compatibility between diamond Ge and 
a phases is of major importance here. We now evaluate 
interface energies within the usual slab approximation, 
ie., from simulations of alternating crystalline layers. 
The surface energy £^s, between two phases 7 and 5, is 
thus defined by 
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where S is the total boundary area and E is the total 
internal energy in the periodic cell, N^ and N^ are the 
number of atoms in the layers of the respective phases, 
E^ and E^ are the formation energies [Eq. (1)]. Ap- 
pendix B details the relation between this Eq. (4) and 
the thermodynamic semi-grand canonical ensemble. At 
equilibrium, the crystalline layers have not exactly the 
same lattice parameters parallel to the interface. Thus, 
to minimize the strain imposed to the layers by the com- 
mon computer cell periodicities, the size of this cell is 
relaxed perpendicularly to the interface during the cal- 
culation of E. Accordingly, the formation energies E^ 
and Ep in Eq. (4) are calculated with the same elastic 
conditions: fixed unit-cell periodicities except for the di- 
rection perpendicular to the interface. The width of the 
layers is increased until we get a convergence for Es bet- 
ter than 2 meV/A^. This condition is achieved when the 
width of each compound is greater than 15 A, i.e., a value 
that is of the same order of magnitude as the nanocolumn 
diameters. 

The main interface orientations, [10 0] and [110], 
have been considered for the Ge/a-Ge2Mn boundary. 
For each of these orientations, all possible atomic ar- 
rangements at the interface have been tested, result- 
ing in 22 different configurations. The calculated in- 
terface energy between a Ge2Mn and diamond Ge is 
found to be about 35 meV/A^ for the lateral bound- 
aries of the column (i.e., [10 0] and [110] orientations. 



the column axes being along [001]). For the inter- 
face perpendicular to the column axis, the interface en- 
ergy is found larger: 63 meV/A^, which is compatible 
with the cylindrical shape of the grains. These ener- 
gies are to be compared with the larger value of Es 
for the interfaces between MusGes and diamond Ge 
(Ref. 24): 80 meV/A^ (Mn5Ge3 [1121] on Ge [1 1 0] and 
); 87 meV/A^ (MnsGes [0 01] on Ge [001]). Only the 
boundary corresponding to MusGea [0 01] on Ge [1 1 1] 
has a similar energy: 35 meV/A^. 

The transformation of the nanocolumns into rounded 
grains after annealing, shows that the nanocolumns are 
out of equilibrium. Above, we have introduced and calcu- 
lated the main terms of energy: chemical energy, elastic 
energy and interface energy. We now want to take into 
account all of these terms and discuss the stability of the 
system depending on the particular distribution of Mn 
atoms. 

First, we consider three types of large scale order: i) 
full homogeneous dilution of Mn atoms in the whole Ge 
diamond matrix; ii) nanocolumns with a high Mn concen- 
tration embedded in a matrix of germanium depleted in 
Mn; and iii) spherical clusters under similar conditions. 
The first case has been assumed in Refs. 8 and 9, with 
a growth at low temperature probably corresponding to 
a very limited atomic diffusion. The second situation 
has been obtained at intermediate temperature^'^^; the 
columnar shape has been linked to a predominance of 
surface diffusion during the growth^^. Finally, the spher- 
ical clusters have been obtained with a high tempera- 
ture growth^^ or after annealing^, and therefore involve 
a three-dimensional (3D) diffusion. 

Second, we compute the total energy for different 
atomic local orders: a) simple defects of the Ge diamond 
crystal; b) Qf-Ge2Mn crystalline grains; or c) MusGea 
crystalline grains. The details of the size and morphol- 
ogy are taken from Ref. 17: average Mn concentration of 
4.1%, 30 000 columns per /im^ (which together fix their 
mean diameter, 2 nm for a Ge2Mn for instance), and 
spherical clusters of 10 nm in diameter. 

Figure 3 presents the energy contributions when the 
different chemical orders are considered. First, it appears 
that the solid solution is ruled out by its very high total 
energy, in both phase-separation cases. The case of full 
dilution has therefore to be limited to lower Mn contents. 
The high surface-to-volume ratio of the nanocolumns 
{S/V = 2 nm~^ here as compared to 0.6 nm~^ for 
a 10-nm sphere) is indeed not high enough to favor a 
columnar solid solution: the low chemical energy of a 
Ge2Mn compensates for its interface energy. On the 
other hand, S/V is yet high enough to make the for- 
mation of MnsGes energetically unfavorable, because of 
its structural incompatibility with diamond- Ge, even if 
its formation energy is almost zero. Indeed, the lateral 
coherence of a Ge2Mn with the Ge crystal, and the re- 
sulting low-interface-energy compensate in this case for 
its higher formation energy. Finally, in the case of spher- 
ical precipitates where the ratio S/V is much smaller. 
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Figure 3. (Color online) Computed energies per atom of a 
Ge-Mn system, embedded in a Ge matrix depending on its 
crystalline structure: Ge-Mn solid solution, a Ge2Mn, and 
?7-Mn5Ge3. The energy model includes chemical, elastic, and 
interface energies. It is applied to three distinct Mn reparti- 
tions, which are met ast able but experimentally observed de- 
pending on the growth conditions: i) Mn dilute in the host 
Ge matrix (DMS), ii) Mn atoms gathered inside nanocolumns, 
and iii) Mn inside spherical precipitates. The assumed char- 
acteristics of the atomic diffusion in these three cases are in- 
dicated. The experimental data describing these Mn repar- 
titions are taken from Devillers et al}^ : global Mn concen- 
tration is 4.1%, column density is 30000 /im~^, diameter of 
spherical precipitates is 10 nm. The zero of the energy cor- 
responds to a phase separation between diamond Ge and 0- 
MnllGeS (the atom reservoirs) with interfaces, free surfaces, 
and elastic constraints that are negligible. 



the interface-energy advantage of a Ge2Mn phase is not 
as crucial, and the MusGes compound is found stable as 
observed experimentally. 

From these observations, we can now propose a coher- 
ent precipitation sequence in the line of the three-stage 
precipitation sequence of AI2CU (diluted -^ 0" -^ 0' ^ 
6). The first stage (DMS -^ 6" ) corresponds to high 
surface- diffusion during the sample growth^ ^ and makes 
possible the precipitation of Mn-rich domains in the sub- 
surface of the sample. Such a structure has not been 
studied in the present paper but is an extension of the 
reported subsurface interstitials^'^^. The second stage 
iO" -^ 9') corresponds to the germination of a coherent 
6' phase from Mn interstitials. This phase has been iden- 
tified as the proposed a phase, which is more likely to be 
stabilized in the form of nanocolumns (Fig. 3). Finally, 
the last stage {0' -^ 9) is the transformation into the 
stable and less coherent MusGes compound thanks to 
3D diffusion. This latter transformation also leads to a 
change of the shape of the precipitates from columns to 
round-shape precipitates. 



V. COMPARISON TO EXPERIMENTAL 
RESULTS 

In previous sections, we have assessed the stability of 
the a phase based on theoretical arguments. We now 
want to analyze our previous experimental characteriza- 
tions of the GeMn nanocolumns^'^^'^^^^ taking into ac- 
count the structure that we propose with the a phase. In 
addition, we first present new results of high-resolution 
electron microscopy. 

a. TEM — Up to this work, the hypothesis made 
for the nanocolumn structure was the solid-solution 
phase in particular because of its strong similarity with 
the germanium matrix when observed by transmission 
electron microscopy (TEM) as shown in Figs. 4a and c. 
Yet, the <0 1> atomic rows of the a Ge2Mn phase 
also form a square lattice compatible with the diamond 
Ge peaks positions in Fig. 4c. The presence of two non- 
equivalent sub-lattices, one for Ge and one for Mn, cor- 
responds to a super-structure which should be visible 
in plane view. However, these sub-lattices generate two 
variants depending on the position of the cell origin with 
respect to the diamond-Ge cell, and these two variants 
are expected to coexist inside each column, because of 
the very low interface energy of the anti-phase bound- 
aries bridging them (we calculate it to be lower than 
10 meV/A^). We have included these two variants to cal- 
culate a TEM image of an unrelaxed a-Ge2Mn nanocol- 
umn in germanium (jEMS software^^), in cross section 
(Fig. 4b) and in plane view (Fig. 4d). The simulated im- 
ages show a striking resemblance between a phase and 
diamond Ge. The global contrast variation between the 
two phases is due to the manganese atoms, the germa- 
nium density being identical in both cases. Our calcula- 
tions reproduce the dark rings around each column, cor- 
responding to the structure discontinuity at the bound- 
ary. 

h. Diffraction — With the same inclusion of the two 
phase variants, we have computed the diffraction pattern 
of a collection of nanocolumns embedded in 100 x 100 unit 
cells domains (see Ref. 27 for details on the column dis- 
tribution and on the method). We find that the inten- 
sity is localized at the corresponding diamond diffraction 
peaks, and that the destructive interferences cancel the 
peaks due to the a-phase sub-lattices. 

c. EXAFS — The Mn local atomic environment 
probed by extended X-ray absorption fine structure (EX- 
AFS) spectroscopy at the Mn K-edge is also compatible 
with the proposed a Ge2Mn phase when considered in the 
form of small nanocolumns (2 nm in diameter) embedded 
in a Ge matrix. To demonstrate this, a different approach 
to the standard EXAFS analysis is adopted. Instead of 
fitting the EXAFS signal with a set of theoretical scatter- 
ing paths coming from Mn simple defects in Ge or clus- 
ters from well-known bulk structures (e.g., MusGea), as 
previously done by our group^^ and others^^, we proceed 
in a qualitative analysis of the simulated EXAFS spectra 
based on molecular dynamics (MD) calculations of large 




Figure 4. Experimental (a and c) and simulated (6 and d) 
TEM images of columns of about 2 nm in diameter, a and h 
show cross section along the <1 1 0> direction, while c and 
d are plane views along the <0 1> direction. Simulations 
were performed on an unrelaxed a-phase nanocolumn using 
the JEMS software^^. 




Figure 5. (Color online) Experimental data (Mn K-edge k^- 
weighted EXAFS) for a sample with nanocolumns (NCs^"^^) 
- Ref. 25 - compared with simulated spectra using the a 
Ge2Mn phase in bulk form (a-bulk**^, top right panel) and 
as nanocolumns embedded in Ge (a-NCs**^, bottom right 
panel). 



relaxed cubic cells containing 96 atoms and 11.315 A of 
side, using the experimental lattice parameter of germa- 
nium'^^. Due to the fact that we are considering complex 
structures with multiple absorbing sites, this approach 
has the advantage to correctly account for both the con- 
figurational and dynamical disorders. In fact, the EX- 
AFS spectra are obtained by averaging the spectra of the 
single Mn atoms in the structure on the last 200 frames of 
the MD calculation^^ and successively averaging over all 
the Mn atoms in the cell. For each configuration, a clus- 
ter of i? = 8 A around a central Mn atom is considered, 
and the simulations are carried out using the multiple 



scattering formalism with the Hedin-Lundqvist exchange 
potential as implemented in the feff8 code'^^. 

In Fig. 5 the simulated EXAFS spectra for the 
a Ge2Mn phase in bulk form (a-bulk**^) and as 
nanocolumns embedded in Ge (a-NCs^*^) are compared 
with the experimental data collected from a sample con- 
taining Mn-rich nanocolumns with a nominal Mn doping 
of 6% - details are given in Ref. 25. Both structures 
reproduce well the frequency of the experimental data, 
due to the value of the main Mn-Ge first shell distance 
of 2.51 A( 2% increase to respect bulk Ge) found by the 
previous calculations for the a phase. There is a con- 
siderable difference in the amplitude of the oscillations 
and their shape in the low-/c region. In fact, the double 
bump at 3 < /c < 4 A~^ that is the fingerprint of the 
nanocolumns^^ and the peak at /c = 5 A~^ are better re- 
produced when moving from the bulk form to embedded 
nano-objects. This is demonstrated by the strong ampli- 
tude reduction when moving from the bulk (a-bulk**^) 
to the ideal nanocolumn (a-NCs^*^). This difference can 
be ascribed to the large number of Mn atoms located 
at the interface with the Ge matrix of which local envi- 
ronment is too highly distorted to generate a coherent 
oscillation. Indeed, for nanocolumns of 2 nm in diam- 
eter, half of the Mn atoms are at the interface and are 
displaced by more than 0.2 A as compared to perfect 
a Ge2Mn. In fact, the first-shell coordination number 
^Mn-i found for a-bulk^*^ is 8.5±0.5 while it is 5.5±0.5 
for a-NCs**^. The strong decrease of nMn-i illustrates 
a new way of interpreting the global coordination num- 
ber of 4.0±0.5 found experimentally, not as a single low 
coordinated site but as a mixture of different Mn envi- 
ronments (in the nanocolumns, at their interfaces and 
diluted in Ge). 

d. XAS, XMCD, and magnetometry — Finally, the 
computed electronic properties of the a phases and their 
ferromagnetism are in line with the metallic character ob- 
served by X-ray absorption spectroscopy (XAS) at Mn 
L edge^^ and the experimental magnetic local moment 
around 1 /ie^'^^ from magnetometry and X-ray magnetic 
circular dichroism (XMCD). In particular, it was shown 
that precise eigenvalues of the initial 2p-states were neces- 
sary to obtain accurate DFT-calculated XAS and XMCD 
spectra^'^. The Curie temperature could not be calcu- 
lated using standard techniques'^ because of the extreme 
variety of environments. However the values of 340 K 
and 1.2 /ie/Mn measured for the high-pressure MnGe4 
compound^ ^ give an idea of the potentialities of these 
phases. 



interstitial sites of a simple-cubic Ge crystal and referred 
to as the a phase. Our results show that the first hypoth- 
esis of a solid solution is not compatible with the experi- 
mental manganese concentration, which is between 15% 
and 40%. On the other hand, the proposed metastable 
a phase is coherent with the diamond structure of the 
matrix with which it has a low interfacial energy. 

The solid solution has a negligible interfacial energy 
with diamond Ge, but a high formation energy. Con- 
versely, the compound MusGea has a negligible forma- 
tion energy, but a high interfacial energy. The growth 
conditions yield columns, which have a high S/V ratio. 
We have shown that with such a geometry, the a-phase 
nanocolumns have a favorable balance between these two 
energy terms. Thus, the geometry imposed by the ma- 
terial growth stabilizes this otherwise metastable com- 
pound. 

The a phase has a coherent interface with the ger- 
manium matrix and has a close connection with a sim- 
ple cubic lattice. Experimentally, it is difficult to dis- 
tinguish it from an hypothetical diamond cubic Ge lat- 
tice with highly concentrated Mn impurities. Although 
our argument in favor of the a phase as a constituent of 
the nanocolumns is based on theoretical results, we have 
shown these results are coherent with the experimental 
data available so far. 

Finally, with the GeMn system as a prototype, we have 
adapted the three-step mechanism of precipitation in the 
AlCu system to rationalize the formation of coherent 
solute-rich precipitates in magnetic semi-conductors, the 
a phase playing the same role as the metastable AI2CU 
6^ phase. Such an analysis would also prevail for oth- 
ers multifunctional semi-conductors where heterogeneous 
materials are used to tailor two or three different prop- 
erties. 
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Appendix A: The manganese point defects in 
diamond germanium 

1. Introduction 



CONCLUSIONS 

We have considered two types of (Ge,Mn) compounds 
as possible candidates for coherent nanocolumn precipi- 
tates in germanium: i) a solid-solution model based on di- 
amond Ge with substitutional and interstitial Mn atoms; 
ii) a new intermetallic compound with the Mn atoms at 



This appendix is the starting point of our analysis of 
the Ge-Mn solid solution done in section II. Indeed, we 
detail here the results regarding the point defects con- 
sisting of dilute Mn atoms in the Ge diamond structure 
and the dimers of such defects. We address in particular 
the issue of their formation energy and their magnetic 
moment. 

We consider the six following defects: 




Figure 6. (Color online) The Mn point defects in diamond 
Ge. Only (met a-) stable configurations are presented. 



• Mns : the substitutional position, which is the mere 
replacement of a Ge atom by a Mn atom; 

• Mut : the tetrahedral interstitial, in which the Mn 
atom is inserted in the bigger hole of the diamond 
structure; 

• Muh : the hexagonal interstitial, for which the Mn 
atom resides at the center of a Ge (non-planar) 
hexagon; 

• Mnv : the split-vacancy, which is the final stable 
configuration of the initial combination of a va- 
cancy and a neighboring substitutional Mn atom; 

• Mud : the "dumbbell" interstitial, for which one 
Mn and one Ge share one lattice site; 

• MnpFCD • the ffcd {fourfold coordinated point de- 
fect^^)^ which is a defect where bonds are rear- 
ranged keeping constant the coordination number 
of each atom. In our case, one of the two central 
atoms composing this defect is a Mn atom. 

We can split these defect into three categories, depending 
on the deviation A7V in the number of atoms compared 
to the pure diamond lattice: i) those with zero variation 
(Mns, Muffcd), those for which one atom is removed 
(Mnv), and those for which one atom is added (Mux, 
Muh, Mud). Among all these defects, two are found 
unstable upon structural relaxation, namely Mud and 
Muffcd- They become Mut and Mns respectively. The 
four others are presented on Fig. 6. 



2. Computational details 

The DFT computation of heterogeneous defects in dia- 
mond Ge is not straightforward for three main reasons: i) 



Defect 


Aiv 


^F (eV) 


Mtot 


C?Mn-Ge 






this work 


Ref. 19" 


Ref. 36" 


(Mb) 


(A) 


Mns 





1.5 


1.6 


2.0 


3 


2.453 (-1.8%) 


Mut 


+1 


2.0 


2.3 


3.0 


5 


2.586 (+3.6%) 


Muh 


+1 


2.3 






5 


2.518 (+5.2%) 


Mnv 


-1 


2.8 






4 


2.761 (-12.1%) 



^ The shown results have been adapted to take into account 
reference energy differences. 

Table I. Main characteristics of the four stable Mn simple de- 
fects in diamond Ge, simulated on Ge's equilibrium lattice 
parameter: formation energy ^f (stabler defect is in bold), 
magnetic moment Mtot and first neighbor distance c/Mn-Ge 
(the deviation from pure diamond first neighbor distance is 
indicated). Aat is the change in atom number in the simula- 
tion cell, as compared to pure diamond. 
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Figure 7. Evolution of the formation energy of the 4 stable Mn 
simple defects versus isotropic strain of the diamond Ge. The 
experimental lattice parameter is indicated by the vertical dot 
dashed line at -1.9%. 



due to the exchange-correlation approximation made for 
the calculation, the diamond-Ge gap is reduced to almost 
zero, and some localized electronic states of the defect 
could interact unrealistically with the diamond Ge con- 
duction band; ii) again due to the exchange-correlation 
approximation, the equilibrium lattice parameter of dia- 
mond Ge is off by almost 2%, which is expected to have 
an impact on the defect characteristics; iii) the compu- 
tation of the formation energy of the defects requires the 
use of reference energies, which is obvious for Ge, but 
not for Mn. We will now briefly detail how we intend to 
overcome these major issues. 

e. Diamond-Ge gap and k -points — To investigate 
the properties of an isolated defect, one must ensure that 
the defect is actually isolated in the simulations. How- 
ever, the standard DFT methods rely on periodic bound- 
ary conditions to be efficient, yielding the simulation of 
not one, but an infinite array of defects. The distance 
between image defects is thus a key parameter that con- 
trols the accuracy of the simulation. Here, we use cubic 
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(eV) 


(Mb) 
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2.00 (-20.0%) 
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2.52 (+0.1%) 


MnsT 


+ 1 
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-0.3 
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2.45 (-2.0%) 


MnD2 


+ 1 


2.4 


-0.1 





2.31 (-7.4%) 


MnFFCD2 


+0 


1.9 


-0.6 





2.34 (-6.5%) 


Mns" 





1.7 


- 


3 
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+4 GPa); +2% (- -3 GPa); and +4% (- -5 GPa). 

g. Formation energy — The formation energy Ep is 
given by a relation similar to Eq. (1), where this time the 
energy is not a quantity by atom but for the defect: 



^ Simulated in a 32 atom box. 

Table II. Main characteristics of five Mn-dimers in diamond 
Ge, simulated on Ge's equilibrium lattice parameter. The 
characteristic of the substitutional defect Mns, computed 
with the same simulation cell, is also given for comparison. 
Aa^ is the change in atom number in the simulation cell, as 
compared to pure diamond. Ep: minimum formation energy 
per Mn (stabler defect is in bold); A^^ : energy difference 
per Mn between ferromagnetic and antiferromagnetic config- 
urations; Mtot: magnetic moment; (iMn-Mn: Mn-Mn distance 
(the deviation from pure diamond first neighbor Ge distance 
is indicated too). 



boxes that contain 64 Ge atom in the diamond structure. 
In this case, the distance between image defects is about 
11.5 A, i.e., large enough to prevent important interac- 
tions through the strain of the matrix. However, as noted 
earlier, electronic interactions may occur via the conduc- 
tion band of the diamond Ge, artificially lowered by the 
exchange-correlation approximation. A way to limit this 
phenomenon is to purposely use a coarse /c-point mesh, 
along with special /c-points'^^, to minimize the disper- 
sive characteristic of the defect states. Sufficient /c-point 
mesh density must however be used to reproduce accu- 
rately the diamond Ge crystal. In the following calcu- 
lations, we use a /c-point mesh that ensures errors lower 
than 5 meV/atom on the energy and 0.05% on the lat- 
tice parameter of the Ge crystal. It corresponds to a set 
of four k points in the 64-atom box (by symmetry, only 
one inequivalent /c-point), and for the two-atom diamond 
unit cell, to a 4 x 4 x 4 Monkhorst-Pack grid^^ with the 
four usual shifts {[^ ^ l]Al ^ ^]A^ I ^]A^ ^ l])- This 
method is particularly efficient regarding the magnetic 
moment, as it gives the correct integer value in fiB unit, 
i.e., the one obtained for bigger boxes, whereas a larger 
mesh, for example, 6 x 6 x 6 in the 64-atom box, gives 
wrong fractional values. 

/. Lattice parameter and strain — The DFT cal- 
culated equilibrium lattice parameter of diamond Ge is 
found larger than the experimental value by 1.9% with 
the usual Perdew-Burke-Ernzerhof (PBE) GGA approx- 
imation done in this article. Hybrid functionals may re- 
produce the lattice parameter with a better agreement 
as shown in Ref. 39. Here, to directly evaluate the im- 
pact of the lattice parameter deviation on the character- 
istics of the defects, we compute them not only with this 
equilibrium lattice parameter, but also with 4 other lat- 
tice parameters shifted by: -4% (~ +10 GPa); -2% (~ 



^F 



^Mn^Mn " ^Ge^Ge 



(Al) 



where £^tot is the DFT total energy of the cell, A^Mn and 
Nqq are the partial numbers of atoms in the cell, and 
MMn and /iGe the respective chemical potentials. This 
formula is quite straightforward when one sees the chem- 
ical potentials as reference energy, as it has been done 
in many previous works. The chemical potential of Ge, 
/iGe, is easy to determine, namely, by considering the en- 
ergy per atom of the pure Ge crystal. However, jUMn has 
been given very different values in the previous works us- 
ing many different references: the isolated Mn atom; the 
Mn atom as substitutional defect^^; the GeMn zinc-blend 
structure; the high temperature Mn 7 phase'^^; or the 
high pressure GeMn B20 phase^^. Some authors do not 
even use references'^^'^^, but just give formation energy 
differences. We have considered none of these methods 
and, as described in Sec. I, we choose diamond Ge and 
MniiGeg as the reference atom reservoirs. 



3. Results 

h. Simple defects — The results obtained for sim- 
ple defects are gathered in Table I. We can see that the 
substitutional Mn (Mns) is the most stable simple de- 
fect in diamond Ge, and exhibits a magnetic moment of 
3 /iB, as already shown by previous works^^. We find the 
value of 1.5 eV for its formation energy, which explains 
the very low solubility of Mn in diamond Ge, and the 
initial slope of £^F(^Mn) in Fig- 2. Magnetic moments 
are in accordance with other theoretical investigations, 
when their simulation boxes are large enough, or /c-point 
mesh coarse enough. Regarding the impact of the choice 
of lattice parameter on our results, the evolution of the 
formation energy versus the isotropic strain applied to 
the diamond is shown in Fig. 7. As expected, we observe 
that the slope of the energy curve versus strain depends 
on the type of defect, namely, if there are more or less 
atoms than in pure diamond Ge. In the substitutional 
case Mns, in which the number of atoms is conserved, 
strain has almost no impact, and the formation energy 
given in Table I is fairly accurate. On the contrary, for 
interstitials (Mux and Muh), a decrease of the lattice pa- 
rameter (to come closer to the experimental lattice pa- 
rameter) increases the formation energy. Oppositely, for 
the split vacancy Mny, the formation energy decreases 
with the lattice parameter. Thus, if Mns appears to be 
definitively the most stable defect, and Mnx stabler than 
Muh, the question of the relative positioning of Mny re- 
mains. 

i. Dimers — To go a step further in the study of 
these defects, we now study the several possible associ- 
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ations of two close Mn atoms, i.e., dimers, in particular 
the following ones: 

• Mn2s * two substitutional Mn at first neighbor sites 
(2 X Mns); 

• Mn2T • two tetrahedral interstitial Mn (2 x Mux); 

• Must • one substitutional Mn and one tetrahedral 
interstitial Mn (Mns+Mnx); 

• MnD2 * the "dumbbell" interstitial, in which both 
atoms of the dimer are Mn, located on both sides 
of a lattice site; 

• MnFFCD2 • the ffcd, in which both atoms com- 
posing the defect are Mn atoms. 

Simulation box are here twice smaller than for simple 
defects, i.e., those of 32 atoms in the pure Ge crys- 
tal. Each dimer was tested in both ferromagnetic and 
antiferromagnetic configurations. Results are gathered 
in Table II. In addition, the dimers Mny+Mnx and 
Mus+Muh are found to be unstable. They relax toward 
Mn2s and Must respectively. 

Two points are worth noticing. First, as already 
proposed^^, the formation energy of Mn defects is re- 
duced when they combine together, in particular from 
two Mns iiito one Must- Second, dimers with the lowest 
energy (below 2 eV/Mn) are antiferromagnetic or ferri- 
magnetic. 



Appendix B: Interface energy 

We detail in this appendix the calculation of the inter- 
face energy as done in section IV. 

In this article, we have considered the thermodynamic 
semi-grand canonical ensemble^'^ to calculate the forma- 
tion energy per atom at T=0 of an alloy 7 of A and B 
atoms as 



the atom reservoirs) ^^: 



EZ = e^ ■ 



XJ^/iA 



X^QJlB 



(Bl) 



where e^ is the DFT energy per atom, jia and iib are the 
chemical potentials of respectively A and B atoms. The 
quantity x\ = 1 — x^^ \s the concentration of A atoms. 

In a similar way, we can compute the total formation 
energy E^ of a given system of Na atoms A and Nb 
atoms B distributed in two phases 7 and 5: 



Ef = E- NaI^a - NbIJLb 



(B2) 



where the internal energy E is the DFT total energy of 
the system. 

Let ii\ and /i^ be the chemical potentials when the 
phases 7 and 5 are at equilibrium (i.e., 7 and 5 represent 



A 
^1% 



(4eT - xle^)/{x% - xl) (B3a) 

(4e^ - xy)l{x'^ - x\). (B3b) 

Note that the chemical compositions of 7 and 5 must be 
different for the denominators in Eq.( B3) to be nonzero. 
We restrict our discussion to this case. By definition, the 
formation energies of 7 and 8 phases are null for (ia = fJ^A 
and /iB = y^B ^^^ then Eq. (Bl) can also be written as 



M 



f^A)^xl{fl% 



E^ - 



Mb) 
Mb) 



(B4a) 
(B4b) 



Now, we want to calculate the interface energy Ei as- 
sociated to the boundary between 7 and 5. Energy Ei 
is defined as the energy difference between the formation 
energy E^ of the combined system and the total energies 
of the separated phases in bulk forms: 



E\ = Ey 



N^El 



■N^E^ 



(B5) 



where N^ and N^ are the respective number of atoms in 
7 and S phases. 

This definition (B5), allows us to compare different 
metastable structures of the interface, for instance, a 
mere juxtaposition of both types of crystal cells or rather 
reconstructions of the interface zone with different local 
atomic concentrations. 

The respective numbers N^ and N^ of atoms in 7 and 
S are related to the numbers Na and Nb by 



Na = xIN^^x^N^ 



Nn = xlN^ ■ 



VA^ 

-x%N' 



(B6a) 
(B6b) 



This system of equations defines the numbers N^ and 
N^ if, and only if, the compositions of 7 and S phases are 
different. This is indeed the case since 7 and 6 correspond 
to the atom reservoirs that fix /i^ and /i^ [Eq. (B3)]. 

From equations (B4) and (B6), we can relate the bulk 
part of the system energy to the deviations of the chem- 
ical potentials: 

N^E^ + N'e'p = Na{^i°a - ma) + NB{fi% " A^b) (B7) 
and replace this bulk energy in Eq.( B5) that becomes 
Ei = E- NaIi\ - NBfil ( V(ma, Mb) ) ■ (B8) 

Therefore, the interface energy Ei of the boundary be- 
tween the two phases is the formation energy ^p of the 
system calculated for chemical potentials at which the 
phases are both at equilibrium so that their bulk forma- 
tion energies vanish. 
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